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ABSTRACT 


Riccati differential equation is one of the most essential tools for modelling many physical situations, such as spring mass 
systems, resistor-capacitor-induction circuits, and chemical reactions among many others. It is applicable in engineering 
and science, and also useful in network synthesis and optimal control. We derived a quarter-step method for the solution of 
RDEs by collocating and interpolating the Laguerre polynomial basis function which does not require starting values 
before they are implemented and they simultaneously generate approximations at different grid points in the interval of 
integration. To show the accuracy and efficiency of our method, five (5) model RDE problems were solved and results 
obtained in terms of the point wise absolute errors shows that the method approximates well with the exact solution. The 


Stability analysis conducted reveals that our method is zero-stable, consistent and convergent. 
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1. INTRODUCTION 


We consider a numerical method for solving general Ricatti differential equation of the form 
y'(t)= a(t)+ db(thy() + c(t y* (1) (1.1) 
Y(ty) = Vo (1.2) 


Where a(t), b(t), c(t) are continuous with c(t)#0 and ¢,, y, are arbitrary constants for y(t) which is an 


unknown function. 


The RDE in (1.1) can also be denoted by the equation below; 


y'(t)= f(t, y) (1.3) 


The conventional linear multistep method associated with equation (1.3) is given in the form 


k k 
DP nei = hy WS nai (1.4) 
i=0 i=0 
Where, @, and WY, are unknown coefficients of the method to be uniquely determined. 
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A differential equation in which the unknown function is a function of two or more independent variable is called 
a Partial Differential Equations (PDEs). Those in which the unknown function is function of only one independent variable 


are called Ordinary Deferential Equations (ODEs). 


Many scholars have worked extensively on the solution of (1.1) in literatures [1-10,12, 14,15,16 and 17]. These 
authors proposed different method ranging from predictor corrector method to block method using different polynomials as 


basis functions, evaluated at some selected points. 


In this paper, we proposed three-step quarter block hybrid methodusing Leguere polynomials as a basis function, 
evaluated at grid and off-grid points to give the needed discrete schemes. These discrete schemes are then combined to 


form the block hybrid method required for implementation. 
2. METHODOLOGY 


In this section, we proposed a computational method for the solution of RDEs of the form (1.1). Themethod doesnot 
require starting values for its implementation and it generates approximations at different grid points in the interval of 
integration. Another advantage of the method is that it is less expensive in terms of function evaluations when compared to 
the LMMs or the Runge-Kutta methods. The major idea in this work is to approximate the exact solution y(x) to (1.1) on 
the partition, Tq) = [@ = Xp < Xy < Xz << Xp < Xp41 <1 < Xy = H] of the quarter-step integration interval [a, b] 


by Laguerre polynomial basis function given by, 


y(x) == wre-*) (2.1) 


n! dx” 


One of the advantages of the Laguerre polynomial (2.1) is that it is orthogonal with respect to the weight function 
w(x) = e* on [0,), Raisinghania (2014). The quarter-step computational method shall be derived using interpolation 
and collocation procedures. Equation (2.1) is interpolated at a grid point; the first derivative of (2.1) is then substituted 
into (1.1) to obtain a differential system which is evaluated at all grid points. Using this technique, in the form of linear 
multistep method, we obtain linear multistep methods which are then put in block form to obtain the new quarter-step 


computational method. 
2.1 Derivation of the Generalized Computational Method 


Let the Laguerre polynomial approximate solution be given by a function of a single variable of the form, 


yx) = att [Ene] = 14+ e- 1) + 2-4 $2) 4 (2.2) 


n! dx” 


Where r and s are the number of collocation and interpolation points respectively. Now, interpolating (2.2) at 


point x,4, = 0 and collocating its first derivatives at points x,4, = 0(1)k — 1, leads to the following system of equations, 
XA=U (2.3) 
Where 


A= [a9 a, a2 ee ul, U= [Ynfrfn+t cy ae i 
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Xn Xn Xn Xm Xe 

0 1 2x, 3x2 ve (N)xN1 
x= 

0 1 2Xn+r Bie —- (N)xner 


Equation (2.3) is then solved for the a;s using Gaussian elimination method. The values of a;-s obtained are then 


substituted into (2.2) to give, after some manipulations, a continuous linear multistep method of the form, 
V(X) = Lixo GX) Inj +h Vino Bi CO frst (2.4) 
Where naj =V(%m+jh) and fri; = Ff neg vay): Again, a@j(x)and £;(x) are unknown function to be 
determined and are expressed as continuous functions of t by writing. 
fee (2.5) 
Thus, evaluating (2.4) at x = Xn44,Xyn42,°t) Xn4y using (2.5), we obtain a set of (r — 1) single schemes which 
can be written in the form of discrete quarter-step computational method as, 
AMY, = EY, + hdf (y,) + hbF (Ym) (2.6) 
Where 
Y= avis Vaee" Vaal 3 Y, = [Yn—(r41) Ya-@r+2) cs 
F(%m) = Pata frte fate] £ On) = Ufn-o-1» fr-c-2) 7 fal” 


A), E,d and b are r X r matrices.Equation (2.6) is referred to as the discrete computational method which gives 


evaluation at different grid points without overlapping. 
2.2 Derivation of Computational Method with Three Partitions 


Let the approximate solution to (1.1) be given by Laguerre polynomial of degree 5, by allowing r +s —1 in 


equation 2.2, that is, 


e~ qn 


y(x) = Vi-o [F< are] = 720 — 1800x + 1200x? — 300x? + 30x4 — x5 (2.7) 
with the first derivative given by, 

y (x) = —1800 + 2400x — 900x? + 120x3 — 5x4 (2.8) 
Substituting (2.8) into (1.1) gives, 


f(x,y) = —1800 + 2400x — 900x? + 120x3 — 5x4 (2.9) 


1 


Now, interpolating (2.7) at point x,,,, = 0 and collocating (2.9) at points x,4,,r = 0 (=) > leads to a system of 


nonlinear equation of the form (2.3), where 


= Tr 
A= [Ap , Qy , Az » a3 14, As] 
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T 
U= Yar fu tns tty Ufa Sofa yt 





720 —1800x, 1200x?2 —300x3 30x —x5 

0 —1800 2400x2 —900x2 120x3 —5x4 
0 —1800 2400x 1 —900x2 , -120x? ; —5x* 
n+ n Té nt+T6 nt+76 
X =|0 —1800 2400x 1 —900x7 , -—120x? , —5x* 
nts nts nts nts 
0 —1800 2400x 3 —900x? 3 -120x? 3; —5x* 3 
nt+T6 nt+T6 n+ nti 
0 —1800 2400x 1 —900x7 , —120x? , —5x* 
n+Z nt+z n+Z ntZ 


Which yields the following matrices 





720 0 0 0 0 0 
0  -—18000 0 0 0 0 Yn 
225 15 5 a fn 
0  -18000 150 — = a 1 
64 512 65536 | [ a2 ree 
0 18000 300 gee cae 2 a f ‘f 
— — =o —— — —— = 1 
16 64 4096 Ay ntg 
F Neches. sabe 2025 405 405 ds fin 
64 512 65536 | 46 f re 
0 18000 600 22 =P 2 ae 
7 ar ano ~ 256 


Solving (2.3) by Gaussian elimination method for the a;'s,j = 0(1)5 and substituting back into the Laguerre 


polynomial basis function gives a quarter-step method of the form, 


1\1 


W(X) = a(x) In + ALI Binay J = 0(=)= (2.10) 


16) 4 
Where the coefficients of y, and ; are give as, 
Yn n+j g& 
A=1 
Bo= — (24576¢5 — 19200t* + 5600t? — 750t? + 45t) 
Bx. =—(—3072¢5 + 2160¢* — 520¢3 + 45t?) 
16 
fi = —=(-6144¢5 + 3840t* — 760t? + 45t?) (2.11) 
8 
Bs = —3072t° + 1680¢* — 280t? + 15¢? 


16 


Bx = ——(—12288t° + 5760¢* — 8803 + 45¢”) 


And t is given by (2.5). Evaluating (2.10) at t = — (=) “gives a discrete quarter-step computational method of 


the form (2.6). 
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10271 24152 __ 53024 Fla __ 44384 























Dah = Tat ee ae gg ge gd ge it gg ne 
Ynta = Int = hh — Bay a 5760 2 aps 11520 nf ae 
Vn yt = Int Geese +B sont a saad am Ta = saan oie 
nak. = Yn a560 salt, + 7 —chf nti + gag = hf nt ~ Foe sa hf nti + S40 hf neh 
nat = Int sa — hf, + AS ONTs +— sf neat or so +— ees 
3. ANALYSIS OF THE METHOD 
3.1 Order of the Method 
The equation (1.4) associated with lineardifference operator Lis defined by: 
Lly(x):h] = Lizo [ajy Gn + jh) — hBiy' %n + jh)] (3.1) 


Where y(x) is an arbitrary test function that is continuously differentiable in the interval [a,b]. Expanding 


y(x, + jh) and y (x, + jh) in Taylor series about x, and collecting like terms in h and y(x) gives: 
L[y(x): h] = coy(x) + cP hy (x) + cP y(xe) + + + CphPy? (x) (3.2) 


Accordingly, following [11] and [14], the differential operator and its associated Linear Multistep Method are said 


to be of order p if: 





Co = Cy = Cp SH Cp = Cp 


1 = 0,64, #0 
The term C;,,, is called error constant and it implies that the local truncation error is given by 
Entk = CpsihP tty? ™* ey) + OChP**) 
3.2 Consistency 


In the spirit of [13], the linear multistep method is said to be Consistent if it has order p => 1. Analysis of our method 


shows that it is consistent since its order p = 5 > 1 (see table 4.6) 
3.3 Zero-Stability 


The linear Multistep Method is said to be zerostable if no root of the first characteristic polynomial has modulus greater 


than one and if every root with modulus is simple. 
The hybrid block method is said to be stable if the root z of the characteristic polynomial p(z), defined by: 
p(R) = det[RA — A’] 
satisfies |R| < 1 and every root with |Z)| = 1 has multiplicity not exceeding two in the limit as n > 0. 


Theorem 1.1 The necessary and sufficient condition for a method to be convergent is for it to be consistent and 


zero stable 
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3.4 Stability of the Computational Method 
3.4.1 Stability of the Computational Method with Three Partitions 


The equation (3.12) when put together formed the block as 


Yn+1 
10 0 0 0y/%n+4 00,0 a) 
010 0 Olly 0000 41°"? 
ntg | _ Yn-2 
00100 00001 
Vn Yn-1 
0 0 0 1 Off is} Jo 0 0 O 1)h, 
000 0 U)¥n4t 00 0 6 1) 


24152 53024 14726 44384 
15 45 45 AS |[fna1 














11 53 19 3238 Whig 
480 5760 11520 5760 || - = 
ap 2 1 1 31 nte 
60 360 1440 360 |If,.3 
9 21 3 51 f 16 
160 640 1280 640 nty 
1 4 7 4 
30 45 360 45 
ih ee 10271 
45 
251 Tn—4 
000 0 —— 
11520]|fn-s 
29 WIf,- 
+n 0 0. 0 —— |)? 
1440 |Ifn-1 
27 II f, 
000 0 Tan 
000 0 y 
360 


normalizing the matrix as h > 0 


p(Z) = [ZA® — A] 


Z0 0 0 07 (0 000 1 2 00 0 -1 

0 Z0 0 0] jo 0 0 0 1] jo z 0 0 -1 

0 0 Z 0 O]-jo 0 0 O 1/=|0 0 ZO -1 

00 0 Z 0] 10 0 0 0 1) jo © 0 Zz -1 

0000 2! loo 0 o tH lo OO Zz -1 
Z*(Z-1) =0 





That is Z, = 1, Z2 = 23 =Z, =Z, =0 
3.5 Convergence 
ALMM is said to be convergent if it is consistent and zero stable. 


From the analysis on zero stability above and in addition to table 4.6, we conclude that our method is both 


consistent and zero stable, hence convergent. 
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4. NUMERICAL EXAMPLES 
The computational method derived shall be applied to some modeled RDEs to test for accuracy and efficiency. 


Problem 4.1 


Consider the following quadratic Riccati differential equation (RDE). 


dy dy 
— = 2y(t) — t7(t 1)— = - 2(t) +1 
(= = 2x) - © +1) 2 = -y@) + 
subject to the initial condition 
y(0) =0 
with theoreticalsolution given as 
2t _ 
j=——— 
ye) ett+] 
Problem 4.2 
Given the RDE 
dy 
—=2y(t)-y*+1 
4 ae 
with initial condition 
y(0) = 0 
and theoretical solution obtained thus, 
1 v2-1 
t)=1+4+v2tanh| V2t + —lo (=) 
y(t) aes eae 


Problem 4.3 
Consider the following nonlinear fractional Riccati differential equation 
D*“y(t) = 1+ 2y(t) — y?(t),0 < On1 
and initial condition 
y(0) = 0 
the theoretical solution for x= 1 was found to be 


V2-1 
V2+1 





1 
y(t) = 1+ V2tanh(v2t + 3 log ); 


Problem 4.4 


Let us consider the problem 
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The theoretical solution 


Problem 4.5 


a =1-u?(x),0<x <1; 


u(0) = 0. 


is 


U(x) = 


e** — J 


e2x* + J’ 


Consider a Riccati differential equation with constant coefficients gien as 


The theoretical solution 


y(t) = y(t) — y(t), 


for this problem is 


et 


~ gee 


y(0) = 0, 


Table 4.1: Numerical Results for Problem 4.1 
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t-value | Theoretical Solution Approximate Abs Error | Time 
0.1000 | 0.42188080599029 | 0.42188080599029 | 1.6653e-16 | 0.0358 
0.2000 | 0.65507140458179 | 0.65507140458179 | 3.3307e-16 | 0.0495 
0.3000 | 0.84956591586548 | 0.84956591586548 | 3.3307e-16 | 0.0631 
0.4000 | 1.01845953531826 | 1.01845953531826 | 4.4409e-16 | 0.0736 
0.5000 | 1.16698055323955 | 1.16698055323955 | 4.4409e-16 | 0.0865 
0.6000 | 1.29820152510130 | 1.29820152510130 | 6.6613e-16 | 0.0950 
0.7000 | 1.41436135877797 | 1.41436135877797 | 6.6613e-16 | 0.1028 
0.8000 | 1.51728483742558 | 1.51728483742558 | 6.6613e-16 | 0.1109 
0.9000 | 1.60853886960997 | 1.60853886960997 | 6.6613e-16 | 0.1207 
1.0000 | 1.68949839159438 | 1.68949839159438 | 6.6613e-16 | 0.1275 
Table 4.2: Numerical Results for Problem 4.2 
t-value | Theoretical Solution Approximate Abs Error | Time 
0.1000 | 0.09966799462496 | 0.09967169539817 | 3.7008e-06 | 0.0418 
0.2000 | 0.19737532022490 | 0.19738963167954 | 1.4311e-05 | 0.0543 
0.3000 | 0.29131261245159 | 0.29134303769860 | 3.0425e-05 | 0.0680 
0.4000 | 0.37994896225523 | 0.37999898473952 | 5.0022e-05 | 0.0803 
0.5000 | 0.46211715726001 | 0.46218801205332 | 7.0855e-05 | 0.0907 
0.6000 | 0.53704956699804 | 0.53714038396459 | 9.0817e-05 | 0.1301 
0.7000 | 0.60436777711716 | 0.60447600054929 | 1.0822e-04 | 0.1404 
0.8000 | 0.66403677026785 | 0.66415872119655 | 1.2195e-04 | 0.1504 
0.9000 | 0.71629787019902 | 0.71642932361034 | 1.3145e-04 | 0.1603 
1.0000 | 0.76159415595577 | 0.76173083991998 | 1.3668e-04 | 0.1697 
Table 4.3: Numerical Results for Problem 4.3 
t-value | Theoretical Solution Approximate Abs Error | Time 
0.1000 | 0.90909090909091 | 0.90902868973486 | 6.2219e-05 | 0.0361 
0.2000 | 0.83333333333333 | 0.83322783539319 | 1.0550e-04 | 0.0459 
0.3000 0.76923076923077 0.76909402754165 | 1.3674e-04 | 0.0576 
0.4000 0.71428571428571 0.71412555961818 | 1.6015e-04 | 0.0646 
0.5000 | 0.66666666666667 | 0.66648827625894 | 1.7839e-04 | 0.0746 
0.6000 0.62500000000000 | 0.62480682274283 | 1.9318e-04 | 0.0853 
0.7000 0.58823529411765 0.58802962114848 | 2.0567e-04 | 0.0950 
0.8000 | 0.55555555555556 | 0.55533888221389 | 2.1667e-04 | 0.1048 
0.9000 | 0.52631578947368 | 0.52608905048065 | 2.2674e-04 | 0.1145 
1.0000 | 0.50000000000000 | 0.49976372539295 | 2.3627e-04 | 0.1223 
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Table 4.4: Numerical Results for Problem 4.4 












































t-value | Theoretical Solution Approximate Abs Error | Time 
0.1000 | 1.12295995501999 | 1.12295995501999 | 5.6155e-13 | 0.0369 
0.2000 | 2.33036366723934 | 2.33036366723934 | 1.1653e-12 | 0.0506 
0.3000 | 3.35929859139219 | 3.35929859139219 | 1.6800e-12 | 0.0642 
0.4000 | 4.07625619989395 | 4.07625619989395 | 2.0384e-12 | 0.0750 
0.5000 | 4.50864023794231 | 4.50864023794231 | 2.2542e-12 | 0.0876 
0.6000 | 4.74705986375187 | 4.74705986375187 | 2.3741e-12 | 0.0972 
0.7000 | 4.87206646548955_| 4.87206646548955 | 2.4363e-12 | 0.1066 
0.8000 | 4.93588015111826 | 4.93588015111826 | 2.4682e-12 | 0.1160 
0.9000 | 4.96801151790818 | 4.96801151790818 | 2.4842e-12 | 0.1268 
1.0000 | 4.98407836223864 | 4.98407836223864 | 2.4922e-12 | 0.1343 











Table 4.5: Numerical Results for Problem 4.5 















































t-value | Theoretical Solution Approximate Abs Error | Time 

0.1000 | -0.09966799462496 | -0.09966799462496 | 2.1917e-010 | 0.0307 
0.2000 | -0.19737532022490 | -0.19737532022490 | 1.6820e-009 | 0.0423 
0.3000 | -0.29131261245159 | -0.29131261245159 | 5.3018e-009 | 0.0549 
0.4000 | -0.37994896225523 | -0.37994896225523 | 1.1443e-008 | 0.0661 
0.5000 | -0.46211715726001 | -0.46211715726001 | 1.9872e-008 | 0.0754 
0.6000 | -0.53704956699804 | -0.53704956699804 | 2.9881e-008 | 0.0852 
0.7000 | -0.60436777711716 | -0.60436777711716 | 4.0496e-008 | 0.0937 
0.8000 | -0.66403677026785 | -0.66403677026785 | 5.0713e-008 | 0.1008 
0.9000 | -0.71629787019902 | -0.71629787019902 | 5.9682e-008 | 0.1064 
1.0000 | -0.76159415595577 | -0.76159415595577 | 6.6810e-008 | 0.1159 











Table 4.6: Order and Error Constants of the New Computational Method 





























Point of Evaluation Cyt Error Constant 
P F 393 
aie 655360 
3 
Xd 5 ee 
iG 2684354560 
ae 5 eee, 
8 1509949440 
3 
x3 5 neers 
tag 2684354560 
zi 5 —— 
31708938240 
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The desirable property of a numerical solution is to behave like the theoretical solution, that is, a good numerical solution 
is one which always converges to its theoretical solution as h > 0. This is an essential property that all numerical methods 
should possess. We are therefore very confident in presenting our results in table 4.1 to table 4.5. In all the examples 
provided, our method tends to converge to its theoretical solution faster, making it one of the favorable methods to be 
considered in solving real life problems. We therefore recommend our method to the scientific world for further 


investigation and application. 
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